Quantum decay rates for driven barrier potentials in the strong friction limit 
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Quantum decay rates for barrier potentials driven by external stochastic and periodic forces in the 
strong damping regime are studied. Based on the recently derived quantum Smoluchowski equation 
[Phys. Rev. Lett. 87, 086802 (2001)] explicit analytical and numerical results are presented for the 
case of the resonant activation phenomenon in a bistable potential and the escape from a metastable 
well with oscillating barrier, respectively. The significant impact of quantum fluctuations is revealed. 
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Introduction- Escape over a high potential barrier 

' driven by thermal noise is of fundamental interest in 

I physics and chemistry And it is meanwhile well 

' understood: after a transient period of time the decay 

[ is governed by a rate constant, which in the simplest 

' classical case is of Arrhenius type provided the decaying 

I system stays in thermal equilibrium far from the bar- 
rier. What happens if due to additional external forces 

] the system is far from thermal equilibrium? This ques- 

■ tion has recently gained much attention for classical bar- 
I rier transport in the strong damping limit (Smoluchowski 

■ limit). Important examples are stochastic resonance [^, 
I resonant activation in biochemical reactions and tunnel 

■ diodes or directed motion in ratchet systems Q. 
] While even the classical physics of such phenomena is 

• in many cases not yet completely explored, much less is 
known about corresponding quantum systems. 

Namely, in contrast to the classical range, tractable 
equations of motions of quantum dissipative systems also 
for strong damping and low temperatures do not exist Q . 

■ The path integral representation provides an exact ex- 
[ pression for the time dependent density matrix, but even 

■ a numerical evaluation is in general prohibitive. Particu- 
I lar progress has been made with the development of Mas- 

■ ter equations |^ and the quasiadiabtic propagator ap- 
proach Q. This way, extensive studies exist for bistable 

■ systems driven by external periodic forces 1^. However, 
' these and related techniques require-either to be valid 

■ or to be practicable-that energy level broadening due 
I to friction remains sufficiently small (depending on the 

• approach) so that for the relevant dynamics the Hilbert 
I space of the bare system can be reduced to a few lowest 

lying eigenstates. Of course, this condition fails in the do- 
main of very strong friction. Recently we showed |^] that 
I exactly in this range crucial simplifications arise: As in 
the classical Smoluchowski limit, momentum equilibrates 
on time scales much faster than any other time scales; ac- 
cordingly, an equation of motion-the so-called quantum 
Smoluchowski equation-can be derived for the position 
probability distribution from the exact path integral re- 
sult. The influence of quantum fluctuations turns out to 
be substantial. Hence, for the first time we are now in 
a position to explore driven barrier escape for arbitrary 



overdamped quantum systems |p^ . 

Here, we focus on two paradigmatic examples, namely, 
a bistable potential with a barrier fluctuating randomly 
in time and a metastable well driven by a periodic 
monochromatic force. Classically, for the first case 
the phenomenon of resonant activation is characteristic, 
while for the second one a substantial rate enhancement 
by driving is observed. We elucidate the significant im- 
pact of quantum fiuctuations on both processes. 

Quantum Smoluchowski-For a systems coupled to a 
heat bath environment the reduced density matrix fol- 
lows from p{t) — Tii,{W{t)} where Trf,{-} denotes the 
trace over bath degrees of freedom and W{t) is the 
time dependent density matrix of the system-|-bath com- 
pound. Within the path integral approach a formally 
exact expression for the position representation p(q, q' , t) 
can be given Now, in Q it was shown that within 
the Quantum-Smoluchowski range (QSR), i.e. 

> nj3, 1/7 and ^7 > fceT, (1) 

the position distribution P{q, t) ~ p{q, q, t) is well de- 
termined by a simple time evolution equation coined 
Quantum Smoluchowski equation (QSE). Here 7 denotes 
the friction constant, luq the ground state frequency in 
a potential V{q), and /3 = inverse temperature. 

This limit is opposite to the classical Smoluchowski range 
where ujq%I3 < 1 and k^T > ?i7. The QSE reads 
P = {l/M-f)dqLq,nP where 

^V' + W" 12 + A:bT9,[1 + \fiV"\ (2) 

with V = dV{q)/dq. Further, 

A = (h/TrM-f) \n{h(3"f/2'K) (3) 

accounts for the dominating impact of quantum fluctu- 
ations. Equivalently, the dynamics of an overdamped 
quantum systems can be seen as a classical Smoluchowski 
dynamics with an effective potential Vqb — V + X/2V" 
and an effective diffusion term Dcs = kBT{l + X/3V"). 
Note that quantum fluctuations in Eq. (|^) are of order 
ln(7)/7 and thus, are much larger than classical finite 
friction corrections which are of order 1/7^. 



Static barriers-We first recall results for the decay rate 
in static barrier potentials. There, adjacent to the bar- 
rier at q = qi, we assume a well region around a min- 
imum at q = qo such that the barrier height VJ, obeys 
Vb S> fceT, huQ. Moreover, we take smooth potentials for 
granted. Then, in |^ it was shown that the escape rate 
out of the well is given within the QSR as 



The substantial rate increase due to quantum fluctua- 
tions well agrees with the exact result (in the limit 
PVb ^ 1). Let us now turn to barrier potentials driven 
by external forces. 

Fluctuating &arriers-Thermally activated diffusion 
over a potential barrier that fluctuates randomly in 
time has evoked much interest recently. In the classi- 
cal Smoluchowski limit it was shown that the interplay 
between relaxation, by thermally activated barrier pas- 
sage, and fluctuation, by correlated external noise, leads 
to a strongly enhanced reaction rate in the resonant ac- 
tivation regime Here, we present the first study to 
this phenomena for corresponding quantum systems. We 
look at the process of dichotomous barrier fluctuations 
with a rate 77 in a symmetric double well and search for 
the ultimate decay rate k{ri) of relaxation to equilibrium. 
Since for this problem an analytical classical theory was 
derived we can adapt the general technique. In par- 
ticular, in the case of high barriers considered here, it 
was shown that rates from the analytical theory are in 
excellent agreement with numerically exact results [ pl| . 

If the potential flips randomly between two surfaces 
V+ and V- at a rate 77, we need two probability densities 
P+{q,t) and P-{q,t) with P+ [P_] being the density to 
find a particle at time t at position q and the potential 
in state V+ [V-^]. Accordingly, the two-dimensional QSE 
reads dtp — with p= (P+,P_) and 




Here, L± — {l/Mj)dqLqm with potentials V± = U ±g. 
The function g describes the barrier modulations and is 
assumed to have the following properties: it is symmetric 
around the barrier top at q — qb and monotone decreas- 
ing away from it; outside some finite range around the 
top it is zero. In particular, = V_ around the well 
minima located at q = igo- Further, its maximum g{qb) 
is small compared to the barrier height Ub but not neces- 
sarily small compared to k^T. The ultimate decay rate 
k(r]) is now defined as the least negative eigenvalue to 



the operator S, 
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FIG. 1 Quantum (solid) and classical (dashed) decay rate vs. 
the barrier fluctuation rate in a double well potential with 
height Ub/hujo — 4 and ujohf3 — 3. Shown are the adiabatic 
approximation Eq. (^) (thick) and the effective potential ap- 
proach Eq. (||) (thin). 

In the adiabatic range of small 77, qualitatively, there 
is no much change to the classical result. After each flip 
the system has enough time to relax to equilibrium inside 
the well regions of the instantaneous potentials. Conse- 
quently, the rate follows as the least negative eigenvalue 
of Sr, when replacing the operators Lqm by the static 
quantum rates — r± for the individual potentials, see 
Eq. (^. Hence, similarly to the classical rate we obtain 

fc(77) = (r+ + r_)/2 + v-[v' + (r+ - r_)V4] . (6) 

As expected k{0) = r+ and k{ri) k^^s = (r+ + r_)/2 
for 77 — » 00. 

More interesting is the region of moderate to fast bar- 
rier flippings, i.e. 77 ^ k{rj). There, the procedure is to 
solve SriP ~ with the boundary condition p{qb) — 
so that due to symmetry p{q) = —p{—q). This way, 
we solve the equation for the equilibrium eigenfunction, 
i.e. with zero eigenvalue, but under boundary conditions 
corresponding to the relaxation eigenfunction with the 
least negative eigenvalue. As long as fc(77) is the smallest 
frequency in the system this approximation is justified. 
Now, in the limit of a high barrier we employ a semiclas- 
sical type of ansatz p = zexp{—(l)/kBT) with an effective 
potential (j) and a prefactor accounting for terms of higher 
order in kBT. To obtain the dominant exponential con- 
tribution we do the following: The two coupled second 
order differential equations corresponding to Eq. (|^) are 
transformed to four coupled equations of first order. This 
linear system is solved by diagonalizing the coefficient 
matrix. The relevant eigenvalue turns out to be — and 
is determined by the solution to the cubic equation 

2A;bT77(1 - XpU")iU' - <f>') - 0'(y| - <P'){VL ~ <P') (7) 

obeying < (t)'U' < {U')^. Here, we introduced effec- 
tive potentials which are related to the bare potentials 
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by y = y + (3A/2)y" - (A/3/2) (F')^ for Y = U,V+, and 
Of course, by formally putting A = in one re- 
gains the known classical result . For the prefactor we 
solve perturbatively the differential equations 5^/? = 0. 
Since this scheme works similarly as in the classical case 
we refer to the literature for further details pl[ |. 

For the rate one first derives from SrjP — —kp that 
m « [P+[qb) + P'_{qb)]/{M^fi![P+ + P-W- Then, in- 
serting the results for P± we gain 



with ljq — U"{qo)/M. The frequency fl{r]) coincides 
with the classical result and thus, its lengthy ex- 
pression is omitted here. However, quantum fluctuations 
strongly affect the dominant exponential factor. In the 
limits of very slow and very fast barrier fluctuations A- 
dependent terms lead basically to a renormalization of 
temperature. For small i] we derive from Eqs. (^) and 
(|) that k{ri) r+. In the region of motional narrowing 
77 — > cx) it is k{ri) — > Tjj where Tu is given by Eq. (^) 
with V replaced by the average potential U. 




FIG. 2 (a) Ratio of the resonant activation rate k^cs to 
fc(0) = F-i- and (b) width of the resonant activation maxi- 
mum vs. A for two different temperatures. Other parameters 
as in Fig. L 

To study the ultimate quantum decay rate for all values 
of r] we evaluated Eq. (||) and Eq. (||) with Eq. (|^) nu- 
merically. Results are displayed in Figs. 1,2. Besides an 
obvious rate increase, the major effects of quantum ffuc- 
tuations are to substantially decrease the relative height 
as well as the width of the resonant activation maximum. 
The shrinking relative height krcs/T+ = (F+ 4- F_)/2F+ 
is a consequence of effectively reduced barrier modula- 
tions g{qb) — > g{qb)[^ — M9"{lb)\]- The decreasing width 
of the plateau range is ascribed for sufficiently large ij to 



the fact that then quantum fluctuations on the left hand 
side in Eq. (0) lead to an effectively enhanced dipping 
rate ij ^ ri[l + Xf3\U" {qb)\] near the barrier top. For in- 
termediate fluctuation rates the cubic equation (|^) gives 
rise to a nonlinear dependence on A leading to an inti- 
mate relation between external barrier fluctuations and 
intrinsic quantum fluctuations. 

Oscillating barriers-A prominent example for a system 
far from equilibrium is a metastable potential driven pe- 
riodically by an external force. Recently, in the classical 
Smoluchowski limit, decay rates were studied for weak 
and moderate to strong driving as well [p^-p^. Here, 
we explore the corresponding quantum problem by cal- 
culating numerically the time dependent decay rate T(t) 
from Eq. for long times. The time dependent poten- 
tial is chosen as V{q,t;a) = Vo{q) + qasm{flt) with a 
static barrier Vo{q) = -(Mtjg/2)g2[i + q/qo]. The static 
barrier top is located at q = 0, the well minimum at 
q = -2go/3, and the static height is Vb/Mujl = 10g§/27. 
For the driven case (a ^ 0) the location of the barrier top 
moves periodically in time with qb{t) ~ (a/Mwf ) sin(r2t). 
Consequently, as in the classical case, the rate T(t) = 

-N{t)/N with N{t) = Jl"^^ P{q,t) also oscillates in 
time. Then, from P — [\/M"f)dqL^^P we derive up 
to negligible corrections 



kBT 
M7 



{1 + xpvi'Mt)]) 



dgP[qb{t),t] 
Nit) 



(9) 



Note that for high barriers and times t <C 1/F(t) one can 
put N{t) « 1 if initially one starts from a normalized 
distribution P{q,0). 




FIG. 3 Time dependent quantum (solid) and classical 
(dashed) rate r(t) for Q/ujt = 2tv (T = 2tv/D.), A = 
a/^hMiol = 0.5, uJbTili = 3, and Vb/Tiujb = 4/3. 

High precision numerical results based on the QSE 
[Eq. (H)] are shown in Fig. 3 for the time dependent rate 
F(<). Apart from a rate enhancement the influence of 
quantum fluctuation is two-fold: first, V{t) tends to be 
more symmetric around its minima and maxima and sec- 
ond, both extrema are slightly shifted to the left. In [ p^ 
Lehmann et al. developed a nice theory for the pure clas- 
sical case (A — 0) based on an asymptotic expansion for 
PVb ^ 1 (for a related theory see also 0). There, one 
finds a weak dependence of the location of the extrema 
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on temperature. Since locally in position space the A 
dependence can be partially incorporated in an effective 
temperature, we attribute the observed shift in the quan- 
tum case to a similar kind of process. However, what is 
not possible is to mimic the numerical data for finite A 
by an effective temperature only. 
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prohibitive. Our intention here is to give for a highly non- 
trivial example a first account about the significant role 
of quantum fluctuations in the QSR. 

To summarize we have studied for two important ex- 
amples quantum transport over driven barrier potentials 
in the strong friction limit. For the case of a fluctuat- 
ing bistable potential we explored the impact of quantum 
fluctuations on the resonant activation phenomenon. In a 
case of a mctastablc potential driven externally by a peri- 
odic force we found a sensitive behavior of the decay rate 
on the relevant parameters. These notables findings may 
stimulate further studies of the quantum Smoluchowski 
equation in physics and chemistry as well. 
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FIG. 4 Arrhenius plot of the time-averaged quantum (solid) 
and classical (dashed) rate T. Other parameters as in Fig. 3. 

The averaged rate 



lim — 



t+T 



dt'm, T = ^ (10) 



as a function either of inverse temperature or driving am- 
plitude is depicted in Figs. 4, 5. Again the rate enhance- 
ment is significant. Interestingly, while on a logarithmic 
scale the classical rate shows a simple linear behavior as 
a function of ujqTi(3 for sufficiently low temperatures (ac- 
cording to the well-known exponential Arrhenius law), 
the quantum rate exhibits a weak nonlinearity. Of course, 
in the high temperature limit quantum corrections be- 
come negligible. As a function of the driving amplitude 
the difference between quantum and classical rate shrinks 
with increasing amplitude. As for strong driving most of 
the escape happens to occur when the barrier is low and 
thus, is much less effective in hindering the transport, the 
effect of a A induced rate enhancement due to a smaller 
barrier height in VeS = V + X/2 V" is diminished. 
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FIG. 5 Time-averaged quantum (solid) and classical (dashed) 
rate as a function of the driving amplitude A = a/ ^ hMujf . 
Other parameters as in Fig. 3. 

What about an analytical theory? We think that in 
principle the classical theory ^3 14 1 can be generalized 
to the quantum detailed analysis, however, still 

needs substantial work. Even then seems an explicit eval- 
uation in closed analytical form for all realistic potentials 
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